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This  report  describes  the  research  performed  under  AFOSR  Grant  No. 
84-0082  during  the  period  February  1,  1986  through  January  31,  1989.  The 
major  objectives  of  the  program  were  (i)  to  determine  the  characteristics  of 
solid  propellant  gas  phase  flames  in  rocket  motors  experiencing  axial 
instabilities  and  (ii)  to  determine  the  validity  of  state  of  the  art  solid 
propellant  response  models.  The  program  was  divided  into  two  tasks  in  order 
to  achieve  these  objectives. 

In  Task  I,  the  response  of  sidewall  stabilized  premixed  flames  to 
longitudinal,  standing  acoustic  waves  (which  simulate  the  oscillations 
encountered  in  unstable  rocket  motors)  was  studied.  A  premixed  flame  was 
chosen  for  this  first  phase  as  it  eliminated  the  need  to  deal  with  difficul¬ 
ties  arising  from  the  presence  of  diffusion  processes  in  the  flame  (these 
were  studied  in  Task  II  of  the  program)  while  providing  a  flame  possessing 
many  important  features  of  actual  solid  propellant  flames.  A  theoretical 
model  of  the  unsteady  behavior  of  such  flames,  based  upon  actual  solid 
propellant  response  models,  was  developed.  Experimental  measurements  of 
velocities  and  heat  release  rates  as  well  as  flow  visualization  studies  were 
undertaken  to  check  the  applicability  of  the  model  and  elucidate+  important 
features  of  the  unsteady  flame  behavior.  It  was  found  that  the  effect  of  the 
flame  driving  or  damping  is  manifested  in  sharp  changes  in  the  normal  veloci¬ 
ty  fluctuation  in  the  flame  region.  This  "pumping"  action  results  in  the 
interchange  of  energy  between  the  core  flow  oscillations  and  the  unsteady 
flame  processes.  It  was  shown,  for  the  first  time,  that  the  flame  behaves  as 
a  combination  of  an  acoustic  monopole  and  an  acoustic  dipole.  Close  agree¬ 
ment  between  theory  and  experiment  was  demonstrated  indicating  that  state  of 
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the  art  models  of  unsteady  solid  propellant  flames  are  capable  of  predicting 
the  flame-acoustic  interactions  during  axial  instabilities.  Furthermore,  it 
was  shown  that  the  acoustic  admittance  at  the  propellant  surface  can  exert 
considerable  influence  on  the  flame  driving  characteristics  and  that  the 
unsteady  heat  transfer  from  the  flame  to  the  propellant  surface  may  provide  a 
feedback  mechanism  between  propellant  burn  rate  and  the  acoustic  oscillations 
which  provides  the  energy  required  to  sustain  axial  instabilities. 

During  Task  II,  the  effect  of  diffusion  limited  processes  on  the  driv¬ 
ing/damping  of  axial  instabilities  in  solid  propellant  rocket  motors  was 
studied.  This  was  accomplished  by  investigating  the  unsteady  behavior  of 
multiple  diffusion  flames  stabilized  on  the  sidewall  of  a  duct.  This  config¬ 
uration  simulated  the  unsteady  burning  of  sandwich  type  solid  propellant 
flames.  It  was  shown  that  such  diffusion  flames  are  responsive  to  axial 
acoustic  oscillations  and  that  the  response  is  frequency  dependent.  A 
theoretical  model  of  the  unsteady  behavior  of  acoustically  excited  diffusion 
flames  was  developed  demonstrating  that  the  effects  of  diffusion  processes 
may  be  incorporated  into  models  of  unsteady  solid  propellant  flame  behavior. 
Finally,  it  was  shown  that  the  behavior  of  the  excited  diffusion  flames  is 
similar  to  the  acoustic  monopole-dipole  combination  behavior  which  was 
exhibited  by  the  investigated  premixed  flames. 
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This  report  summarizes  the  research  performed  during  the  period  February 
1,  1986  through  January  31,  1989  under  AFOSR  Grant  No.  AF0SR-84-0082  which 
is  entitled  "Investigation  of  the  Flame-Acoustic  Wave  Interaction  During 
Axial  Solid  Rocket  Instabilities".  The  overall  objectives  of  this  research 
were  (i)  the  determination  of  the  characteristics  of  solid  propellant  gas 
phase  flames  in  rocket  motors  experiencing  axial  instabilities  and  the 
contributions  of  these  flames  to  the  driving  of  the  instabilities,  and  (ii) 
the  determination  of  the  validity  of  state  of  the  art  solid  propellant 
response  models.  The  problem  is  of  iftuch  practical  interest  because  it  is  well 
known  that  the  interaction  of  solid  propellant  flames  with  the  rocket  flow 
oscillations  is  the  mechanism  responsible  for  driving  the  instability. 
Consequently,  developing  an  understanding  of  the  fundamental  processes  which 
control  this  interaction  may  lead  to  the  development  of  practical  means  for 
reducing  the  driving  provided  by  oscillatory  solid  propellant  flames. 

The  characteristics  of  solid  propellant  flames  depend  upon  the  interac¬ 
tion  between  processes  which  occur  in  the  gas  and  solid  phases.  During 
burning,  heating  from  the  gas  phase  leads  to  propellant  surface  pyrolysis. 
The  gases  leaving  the  surface  continue  reacting,  in  diffusion  and  premixed 
flames,  ’  as  they  move  away  from  the  pyrolyzing  surface.  Some  of  the  re¬ 
leased  thermal  energy  is  transported  back  to  the  propellant  surface  to 
continue  the  burning  process.  The  remainder  of  the  released  energy  causes  an 
increase  in  the  temperature  of  the  gases  from  an  estimated  900°K  at  the 
propellant  surface  to,  approximately,  3300°K  at  the  outer  edge  of  the  flame. 
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In  addition,  momentum  exchange  with  the  axially  moving  core  flow  turns  the 
gases  leaving  the  surface  in  the  direction  of  the  combustor  core  flow. 
Consequently,  the  propellant  burn  rate  depends  upon  gas  and  solid  phase 
chemical  kinetics,  gas  phase  diffusion  and  mixing  processes,  gas  and  solid 
phase  transport  processes  and  so  on. 

It  is  evident,  therefore,  that  the  complexity  of  the  solid  propellant 

flame  presents  the  investigator  with  a  myriad  of  challenging  problems  which 

are  compounded  when  the  multidimensionality  of  the  flow  and  the  unsteady 

effects  are  taken  into  consideration.  Moreover,  even  if  theoretical  models  of 
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these  processes  are  developed,  they  are  not  amenable  to  direct  experimen¬ 
tal  verification.  The  problem  lies  in  the  extremely  small  dimensions  of  the 
solid  propellant  combustion  zone.  For  example,  the  thickness  of  the  solid 
phase  thermal  wave  is  of  the  order  of  15  microns  and  that  of  the  gas  phase 
reaction  zone  is  of  the  order  of  50  microns.1  Probing  such  thin  zones  with 
physical  probes,  or  even  available  optical  probes  is  not  feasible.  These 
difficulties  are  compounded  by  the  short  residence  times  of  the  material 
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within  the  combustion  zone  which  are  of  the  order  of  10  -10  seconds.  In 

addition,  the  smoky  nature  of  the  flame  will  hinder  optical  access  to  the 

propellant  surface  where  combustion  occurs. 

To  deal  with  this  problem,  this  research  program  was  divided  into  two 
8  9 

tasks  ’  '  In  Task  I,  the  response  of  premixed  flames  stabilized  on  the 
sidewall  of  duct  to  the  excitation  of  a  longitudinal,  standing  acoustic  waves 
(which  simulate  the  behavior  of  the  oscillations  in  solid  propellant  rockets 
experiencing  axial  instabilities)  was  studied.  A  premixed  flame  was  chosen 
for  Task  I  of  this  study  because  it  eliminated  the  need  to  deal  with  diffi¬ 
culties  arising  from  the  presence  of  diffusive  processes  in  the  flame  (these 
were  studied  in  the  second  phase  of  the  program)  while  providing  the 
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investigator  with  a  flame  possessing  many  of  the  important  features  of  actual 
solid  propellant  flames.  For  example,  it  possesses  the  sharp  temperature  rise 
between  the  propellant  surface  and  the  flame  edge,1  it  interacts  with  temper¬ 
ature  and  velocity  acoustic  boundary  layers  which  also  exist  next  to  burning 
solid  propellant  surfaces  and  it  provides  a  situation  in  which  an  oscillato¬ 
ry,  multidimensional  flame  region  interacts  with  one  dimensional  core  flow 
oscillations.  In  addition,  the  thicknesses  of  the  investigated  premixed 
flames  could  be  controlled  experimentally  which  provided  one  with  an  opportu¬ 
nity  to  perform  the  needed  measurements. 

The  experimental  set  up  used  during  Task  I  of  this  research  program  is 
shown  schematically  in  Fig.  1.  The  investigated  premixed,  flat  flame  is 
stabilized  above  a  porous  plate  burner  on  the  lower  apparatus  wall  which 
simulates  a  solid  propellant  surface.  Acoustic  drivers  at  one  end  of  the 
apparatus  are  used  to  establish  a  standing  acoustic  wave  of  desired  ampli¬ 
tude,  frequency  and  pressure  and  velocity  relationship  in  the  vicinity  of  the 
flame.  The  behavior  of  the  flame  was  investigated  both  theoretically  and 
experimentally  in  order  to  determine  the  nature  of  its  interaction  with  the 
adjacent  oscillatory  flow. 

The  investigation  specifically  addressed  the  following  two  questions: 

1.  Are  state  of  the  art  models  of  unsteady  solid  propellant  flames  capable 
of  predicting  the  characteristics  of  the  developed  premixed  flame  under 
conditions  simulating  those  encountered  in  unstable  rocket  motors? 

2.  What  features  of  the  flame(e.g. .maximum  flame  temperature,  heat  transfer 
to  the  propellant  surface,  spatial  and  temporal  distributions  of  the 
flame  heat  release  rate  and  so  on)  exert  the  greatest  influence  upon  the 
flame  dri ving/damping  of  the  core  flow  axial  oscillations? 
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Experimental  Set-Up  Schematic 


These  questions  were  answered  by  a  combination  of  theory  and  experiment. 
A  theoretical  model  of  the  premixed  flame  experiencing  axial  acoustic  oscil¬ 
lations  was  developed1^  on  the  basis  of  current,  state  of  the  art  unsteady 
solid  propellant  flame  analyses.  The  predictions  of  this  model  were  compared 
with  measurements  of  the  spatial  distributions  of  normal  velocity  and  heat 
release  rate  and  their  phase  relationships  with  the  acoustic  pressure  oscil¬ 
lations.  As  will  be  discussed  in  the  next  section,  these  quantities  are 
directly  indicative  of  the  acoustic  driving/damping  characteristics  of  the 
flame.  It  was  found  that  the  model  predictions  were  in  excellent  agreement 
with  the  experimentally  measured  data.  In  particular,  the  model  was  able  to 
predict  both  the  frequency  dependence  as  well  as  the  magnitude  of  the  flame 
driving.  An  important  new  ’finding  was  that  acoustic  driving  occurs  by  a 
combination  of  a  monopole  and  dipole  type  acoustic  source  behavior  of  the 
flame.  This  was  corroborated  by  both  theory  and  experiment.  These  findings 
and  their  implications  are  described  in  the  next  section.  Obviously,  these 
studies  greatly  enhance  the  confidence  in  the  validity  of  the  developed  model 
and,  therefore,  in  related  models  of  unsteady  solid  propellant  flames. 

Although  this  agreement  is  most  encouraging,  it  should  be  realized  that 
diffusion  processes  which  are  an  important  facet  of  solid  propellant  combus¬ 
tion  were  not  accounted  for  in  the  first  phase  of  this  program.  These  pro¬ 
cesses  occur,  in  part,  due  to  propellant  heterogeneity1  and  they  involve  the 
diffusion  of  oxidizer  and  fuel  towards  each  other.  These  diffusion  processes 
are  important  as  they  often  control  the  overall  reaction  rates.11  Therefore, 
the  studies  initiated  in  the  first  phase  of  this  research  program  were 
extended  to  investigate  the  behavior  of  carefully  controlled  diffusion  flames 
stabilized  in  longitudinal  acoustic  fields. 
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Task  II  of  the  research  also  involved  a  close  interaction  between 

theoretical  and  experimental  efforts.  The  experimental  set  up  was  modified 

by  replacing  the  side  wall  porous  plug  burner(  see  Fig.  1)  by  a  row  of 

parallel,  alternating  fuel  and  oxidizer  inlet  ports  simulating,  for  example, 

a  sandwich  type  of  propellant  in  which  the  oxidizer  and  binder  portions 

* 

alternate  (see  Fig.  2)  .  Longitudinal  sound  fields  excited  by  the  acoustic 
drivers  were  imposed  upon  the  established  diffusion  flames  to  simulate  the 
unsteady  environment  in  unstable  rocket  motors.  As  with  the  premixed  flame 
set  up,  this  configuration  possessed  many  of  the  features  of  actual  solid 
propellant  flames.  It  simulated  the  sharp  temperature  rise  between  the  solid 
propellant  surface  and  the  flame  edge,1  it  interacted  with  the  temperature 
and  velocity  acoustic  boundary  layers  which  also  exist  next  to  the  burning 
solid  propellant  surface  and  so  on.  In  addition,  the  important  effects  of 
diffusion  processes  in  the  flame  region  were  included  in  the  investigation. 

The  diffusion  flame  studies  were  conducted  with  the  objective  of  answer¬ 
ing  the  following  questions: 

1.  How  does  the  presence  of  diffusion  processes  in  the  flame  region  affect 
the  acoustic  driving/damping  characteristics  of  the  flame? 

2.  What  features  of  the  flame  (e.g.,  the  fuel /air  ratio,  flow  rates,  flame 
temperatures  and  so  on)  exert  the  greatest  influence  upon  the  flame 
driving/damping  characteristics? 

How  may  the  effect  of  diffusion  processes  be  incorporated  in  theoretical 
models  of  unsteady  solid  propellant  flames? 


*Such  sandwiches"  of  actual  solid  propellants  have  been  utilized  in  the  past 
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to  study  solid  propellant  flames. 
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GAS  PHASE  DIFFUSION  FIAME 
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OXIDIZER 

(a) 

DIFFUSION  FIAME 


Fig.  2.  (a)  Sandwich  Propellant  with  Established  Gas  Phase 
Diffusion  Flames  and  (b)  Proposed  Simulation  by 
Gaseous  Diffusion  Flames. 
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In  the  next  section,  the  results  obtained  under  both  tasks  of  this 
program  are  described.  The  significance  of  these  results  is  evaluated  and 
important  findings  are  summarized.  Finally,  Appendix  A  lists  the  personnel 
involved,  the  graduate  degrees  awarded,  and  the  professional  interactions 
made  possible  by  this  program. 

RESEARCH  ACCOMPLISHMENTS 


This  section  describes  the  results  obtained  during  the  two  tasks  of  the 
research  program. 


Task  I: 

Interaction  of  Premixed  Flames  with  Axial  Acoustic  Fields 


The  developed  experimental  set  up  (see  Figs.  1  and  3)  allows  the  pre¬ 
mixed  flame  to  be  stabilized  sufficiently  away  from  the  wall  so  that  the 
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required  measurements  may  be  performed'  It  is  described  in  detail  in 
Ref.  18.  Briefly,  the  set  up  consists  of  a  long  rectangular  tube  having  a 
flat  flame  burner  on  one  of  its  side  walls,  an  axially  movable  injecter  plate 
at  the  inlet  end  and  two  acoustic  drivers  at  the  exhaust  end. 

Ouring  an  experiment,  a  combustible  mixture  of  propane  and  air  is  fed 
into  the  side  wall  burner  and  a  flat  flame  is  stabilized  a  short  distance 
(i.e.,  5-15  mm)  away  from  the  burner  surface  (see  also  Fig.  4).  Next,  the 
acoustic  drivers  are  turned  on  to  excite  a  standing  longitudinal  acoustic 
wave  of  desired  frequency  and  amplitude  in  the  tube.  The  position  of  the 
flame  relative  to  the  standing  acoustic  wave  (i.e.,  next  to  a  pressure  node 
or  pressure  antinode)  can  be  varied  by  moving  the  injector  plate  axially. 
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As  the  results  obtained  during  this  investigation  are  best  described  by 
comparing  experimental  and  theoretical  results,  the  theoretical  model  (de¬ 
scribed  in  detail  in  Ref.  10)  will  be  briefly  outlined.  The  flow  variables 
including  temperature,  pressure,  velocity,  fuel  mass  fraction  and  so  on  are 
expressed  as  sums  of  steady  and  unsteady  components.  A  set  of  nonlinear 
differential  equations  is  obtained  for  the  steady  state  components  and  a  set 
of  linear  differential  equations  for  the  unsteady  components.  The  analysis 
for  the  unsteady  components  is  two  dimensional  and  depends  both  upon  the 
axial  location  (x)  and  the  normal  distance  (y)  from  the  burner  surface.  The 
equations  for  the  unsteady  components  contain  coefficients  which  depend  upon 
the  steady  state  solutions  and  are  subject  to  boundary  conditions  at  the 
burner  surface  (y  =  0)  and  at  the  downstream  edge  of  the  flame  (y  =  y^). 
Thus,  the  steady  state  solutions  need  to  be  determined  either  theoretically 
or  experimentally  or  by  some  combination  of  the  two.  In  addition,  some  of 
the  boundary  conditions,  in  particular  the  value  of  the  normal  component  of 
the  velocity  fluctuation,  v',  at  the  burner  surface  need  to  be  input  from 
experimental  measurements. 

The  acoustic  driving/damping  characteristics  of  the  investigated  flames 
may  be  determined  in  two  ways.  First,  the  acoustic  energy  input  (or  extract¬ 
ed)  by  the  flame  from  the  core  flow  oscillations  is  given  by  the  time  average 
of  the  product  p'v1  of  the  acoustic  pressure  p'  and  the  normal  component  of 
the  velocity  fluctuation  v‘.  In  this  connection,  it  should  be  pointed  out 
that  a  normal  component  of  the  fluctuating  velocity  v'  is  formed  in  the  flame 
region  as  a  result  of  the  periodic  reaction  rate  (and  expansion)  of  the 
flame.  If  the  product  p'v'  is  positive,  then  energy  is  fed  into  the  acoustic 
oscillations  and  vice  versa.  If  all  phases  are  referred  to  that  of  the 
acoustic  pressure  (as  was  done  in  the  experiments)  then  the  pressure 
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oscillations  may  be  taken  to  be  purely  real.  In  such  a  case,  the  energy 

input  of  the  flame  is  given  by  <  p'  Real  (v')>  which  is  the  time  average  of 

the  product  of  p'  and  the  real  part  of  v1  generated  by  the  flame.  This 

yields  a  quantitative  measure  of  the  flame  driving/damping  characteristics. 

Secondly,  a  qualitative  measure  of  the  driving  characteristics  of  the  flame 
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may  be  obtained  using  Rayleigh's  criterion  which  states  that  if  the  un¬ 
steady  heat  release  from  the  flame  is  in  phase  with  the  local  pressure 
oscillations,  then  acoustic  driving  by  the  flame  results  (and  vice  versa)> 
Both  of  these  criteria  were  utilized  in  the  conducted  studies. 


Comparisons  Between  Theory  and  Experiment 
The  model  for  the  unsteady  flame  behavior  requires  the  following  input: 
(i)  steady  state  spatial  distributions  of  the  temperature,  velocities  and 
heat  release  rate. 

(ii)  the  value  of  the  normal  velocity  fluctuation  v1  at  the  burner  surface 
y  =  0. 

In  accordance  with  the  statements  made  earlier  about  the  driving/damping 
characteristics  of  the  flame,  attention  will  be  focussed  herein  on  the 
following  predictions  of  the  model: 

(i)  the  distribution  of  the  real  part  of  the  normal  velocity  fluctuation, 
Real  (v‘ )  and 

(ii)  the  distribution  of  the  heat  release  rate  fluctuation,  q'. 

The  required  steady  state  temperature  and  velocity  distributions  were 

obtained  experimentally.  The  steady  state  temperature  distribution  was 
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determined  us^ng  the  inclined  slit  method  ’  while  the  velocities  were 
measured  using  LDV.  The  steady  state  reaction  rate  (w)  is  related  to  the 


steady  state  temperature,  T,  by  the  following  Arrhenius  type  relation 

(1  -  f )J  ,,, 

W  =  A- — r~  e'E,T 
T 
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where  E  is  the  (normalized)  activation  energy  and  A  is  the  steric  factor.  As 
these  quantities  depend  upon  the  overall  reaction  characteristics  of  the 
flame  and  are  not  readily  available,  the  theoretical  model  was  used,  in  the 
inverse  mode,  to  determine  these  quantities;  that  is,  instead  of  solving  for 
T  using  knowledge  of  A  and  E,  A  and  E  were  determined  from  the  steady  energy 
equation  using  the  experimentally  determined  distribution  of  T. 

A  typical  measured  temperature  profile  and  a  corresponding  steady  state 
solution  obtained  theoretically  are  shown  in  Fig.  5.  The  temperature  is 
plotted  as  a  function  of  the  normal  distance  y  from  the  burner  surface.  Two 
factors  are  of  interest;  the  flame  standoff  distance  is  of  the  order  of  1  cm 
thus  enabling  experimental  probing  and  the  steady  state  temperature  gradient 
is  very  small  in  the  neighborhood  of  y  =  0.  This  means  that  heat  transfer  (by 
conduction)  to  the  burner  surface  is  not  an  important  factor  in  the  experi¬ 
ments.  It  should  be  noted,  however,  that  this  is  not  true  for  actual  solid 
propellant  flames  for  which  the  stand  off  distance  is  of  the  order  of  50 
microns.  The  implication  of  this  will  be  discussed  shortly.  One  should  also 
note  the  sharp  temperature  rise  in  the  flame  or  reaction  zone.  The  corre¬ 
sponding  steady  state  reaction  rate  profile  is  shown  in  Fig.  6. 

The  value  of  v'  at  y  =  0  was  obtained  in  terms  of  the  acoustic  admit¬ 
tance  R  (i.e.,  R  =  v'/p'  at  y  =  0)  of  the  side  wall  burner  surface.  This 

admittance  was  obtained  experimentally  using  the  impedance  tube  technique. 

20  21 

The  measured  admittance  ’  is  plotted  as  a  function  of  frequency  in  Fig.  7 
and  is  seen  to  depend  strongly  upon  it. 

To  determine  the  flame  driving  characteristics,  consider  first  the 
comparisons  of  predicted  Real  (v')  distributions  with  values  obtained  experi¬ 
mentally  using  an  LDV  system.  As  noted  earlier,  this  gives  a  quantitative 
measure  of  the  driving/damping  by  the  flame.  In  Fig.  8,  Real  ( v ' )  is  plotted 
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Typical  Temperature  Distribution  as  a  Function  of  the  Normal 
Distance,  y  (Non-Dimensional),  from  the  Burner  Surface. 


as  a  function  of  y  at  a  frequency  of  200  Hz  (the  driving  frequency).  At  y  = 
0,  it  is  given  by  the  sidewall  admittance  and  is  negative  indicating  that  the 
sidewall  acts  as  an  acoustic  damper.  It  varies  slowly  between  the  burner 
surface  (y  =  0)  and  the  flame  region  where  the  strongest  temperature  rise 
occurs  (compare  with  Fig. 5).  In  this  flame  region  Real  ( v 1 )  becomes  sharply 
less  negative(or  equivalently  more  positive  and  more  in  phase  with  the 
pressure  oscillations).  This  means  that  at  200  Hz,  the  flame  inputs  energy 
into  the  core  flow  oscillations.  Note  also  that  the  prediction  of  the  model 
(the  solid  line)  agrees  extremely  well  with  the  measurements.  It  is  also 
important  to  note  that  although  the  flame  inputs  energy  into  the  acoustic 
field  (by  decreasing  the  "negativeness"  of  Real  (v1)),  it  is  not  able  to 
overcome  the  strong  damping  effect  of  the  sidewall  burner  surface.  This 
implies  that  in  a  solid  rocket  motor  situation,  the  admittance  at  the  propel¬ 
lant  surface  is  an  important  parameter  relating  to  axial  instabilities. 

Consider  next  a  case  where  the  driving  was  at  400  Hz  (Fig.  9).  In  this 
case  Real  (v1)  becomes  even  more  negative  in  the  flame  region  which  indicates 
damping  by  the  flame.  What  is  important  here  is  that  the  model  agrees  with 
the  measurements  and  therefore  demonstrates  its  capability  to  distinguish 
between  situations  where  the  flame  drives  or  damps. 

From  the  Rayleigh  criterion  point  of  view,  the  phase  of  the  unsteady 

heat  release  needs  to  be  compared  with  that  of  the  pressure.  The  heat  release 

fluctuations  have  a  component  in  phase  with  the  pressure  oscillations  if  the 

phase  difference  between  the  two  is  less  than  +  90°.  The  predictions  of  the 

model  are  compared  with  experimentally  measured  values  of  the  phase  in  Fig. 

10.  These  were  obtained  by  measurements  of  overall  CC  radiation  emission  from 
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the  flame.  As  noted  in  earlier  reports,  ’  ’  these  emissions  are  indicative 
of  the  heat  release  rates  from  the  flame.  Note  that  at  200  Hz,  p'  and  q'  are 
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18. 


Steady  State  Reaction  Rate  Profile 


Function  of  Frequency 
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Measured  and  Predicted  Values  of  Real  (v')  as  a  Function 
of  the  Normal  Distance  from  the  Burner  Surface  at  200  Hz 
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Measured  and  Predicted  Values  of  Real  (v  )  as  a  Function 
of  the  Normal  Distance  from  the  Burner  Surface  at  400  Hz 


MEASURED 


Fig,  10.  Measured  and  Predicted  Phase  of  the  Unsteady  Heat  Release 

(Radiation)  with  Respect  to  the  Imposed  Pressure  Oscillations 


in  phase  (indicating  flame  driving  of  the  acoustics)  while  at  400  Hz  they  are 
out  of  phase  (indicating  flame  damping  of  the  acoustics).  These  predictions 
are  in  agreement  with  the  trends  exhibited  by  Real  (v1)  in  Figures  8  and  9. 
Also,  the  theoretically  predicted  trends  are  in  good  agreement  with  the 
experimental  measurements. 

These  comparisons  were  made  with  the  flat  flame  stabilized  at  a  pressure 

antinode  of  the  excited  standing  waves.  Similar  results  were  obtained  with 

the  flame  stabilized  at  other  locations  of  the  standing  wave  as  long  as  this 

location  was  not  at  a  pressure  node.  At  a  pressure  node,  depending  upon  the 

excitation  levels,  the  flame  surface  would  get  distorted  by  the  appearance  of 

spikes  or  wavelets  and  in  extreme  cases  the  flame  would  break  up  and  be 

extinguished.  This  feature  has  been  described  in  earlier  reports  in  detail. 
Only  at  very  low  levels  of  excitation  could  a  stable  flame  be  maintained. 

However,  under  these  conditions  measurements  could  not  be  made  with  any 
degree  of  precision.  Consequently,  efforts  have  been  concentrated  on  under¬ 
standing  the  flame  behavior  away  from  a  pressure  node. 

Important  discoveries  have  been  made  by  considering  the  time  dependent 
motion  of  the  flame.  By  means  of  high  speed  cinematography  (described  in 
References  9  and  20),  it  was  determined  that  under  the  influence  of  a  sound 
field  flames  located  away  from  a  pressure  node  exhibit  an  up  and  down  motion 
relative  to  their  steady  state  location  at  the  frequency  of  the  excited  wave. 
At  the  phase  of  maximum  pressure,  the  flame  would  be  located  closest  to  the 
sidewall  and  at  the  phase  of  minimum  pressure  it  would  be  located  farthest 
from  the  sidewall.  A  prediction  of  the  flame  motion  obtained  using  the 
developed  theoretical  model  is  presented  in  Fig.  11.  It  plots  the  instantane¬ 
ous  heat  release  rate  at  different  times  during  a  cycle  of  oscillation  as  a 
function  of  the  normal  distance  y  from  the  sidewall.  At  any  given  time,  the 


11 


zed ) 
eat 

w 

•  r-4  U 

3 

C3 

O 

CT3  X 

03 

E 

c 

^  03 

03 

o  ^ 

VJ 

3£  03 

C 

'w'  0) 

03 

f— l 

4J 

w  qj 

W 

a  o- 

C 

u 

t-c 

& « 

03 

03 

j£ 

«d 

4-J 

03  £ 

(A 

X  3 

a> 

£ 

c 

4-1  -H 

•r4 

C  X 

£ 

J-i 

03 

c 

4J 

03  03 

03 

o-x 

<u  u 

Q 

Q 

•» 

4-1 

4-1 

• 

03  O 

c 

£ 

** 

o 

•r4  C 

03 

H  0 

E 

4J 

•rH 

•H 

03 

•a  ^ 

M 

a 

<V  <0 
*J  o 

C 

3  3 

03 

> 

<u 

X) 

•r4 

I 

03  0> 

o 

£6 

to 

• 

H 

H 

• 

00 

t4 

C*4 

I — 

<  DU 

S  2 


LJ  O 
i-  on 
01  o 


LlI 

O 

<  ft 

Ll  O 


b  ( t  M  +  m) 


location  of  maximum  heat  release  rate  may  be  identified  as  the  flame  loca¬ 
tion.  As  the  pressure  was  taken  to  vary  as  cos  2Trt/T,  where  T  is  the  period 
of  excitation  and  t  is  the  time,  t  =  0  represents  the  phase  of  maximum 
pressure.  As  is  seen  from  the  figure,  at  this  time  the  flame  is  closest  to 
the  side  wall.  One  half  period  later  it  is  farthest  away  from  the  wall.  Thus, 
in  this  respect  too,  the  predictions  of  the  model  are  in  complete  agreement 
with  the  experimental  observations. 

Much  more  can  be  learned,  however,  from  a  closer  inspection  of  Fig.  11. 
Note  that  when  the  flame  is  closest  to  the  wall,  the  instantaneous  heat 
release  rate  is  also  the  highest  during  the  cycle  and  when  it  is  farthest  the 
instantaneous  heat  release  rate  is  the  lowest.  Similar  results  were  obtained 
experimentally  also  (see  Ref.  9)  in  complete  agreement  with  the  theoretical 
model.  It  was  noted  earlier  that  in  the  experiments,  heat  transfer  from  the 
flame  to  the  sidewall  was  not  an  important  factor  due  to  the  flame  standoff 
distance  being  of  the  order  of  1  cm  which  resulted  in  a  very  low  steady  state 
temperature  gradient  at  the  wall.  In  an  actual  solid  propellant  flame  stand¬ 
off  distances  are,  however,  of  the  order  of  50-100  microns  so  that  heat 
transfer  to  the  surface  of  the  propellant  is  an  important  issue.  If  the  heat 
release  rate  is  highest  when  the  flame  approaches  the  propellant  surface,  as 
is  the  case  with  the  premixed  flame,  a  mechanism  is  available  to  sustain  a 
non  steady  burn  rate  of  the  propellant  in  tandem  with  the  pressure  oscilla¬ 
tions  which  may  lead  to  strong  instabilities. 

The  Flame  Driving  Mechanism 

Up  to  this  point  the  discussion  has  centered  upon  the  excellent  agree¬ 
ment  between  the  developed  model  and  the  experimental  data  which  indicates 
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that  state  of  the  art  models  of  unsteady  solid  propellant  flames  may  indeed 

be  capable  of  capturing  the  salient  features  of  the  flame  behavior.  Now 

attention  is  focused  on  the  driving  characteristics  (i.e.,  consideration  of 

the  flame  as  an  acoustic  source)  of  the  investigated  premixed  flames.  It  is 

22 

well  known  that  the  most  fundamental  acoustic  source  is  a  monopole  which 
may  be  envisioned  as  a  periodically  expanding  and  contracting  balloon.  Other 
acoustic  sources  are  obtained  by  suitable  combination  of  monopole  sources; 
that  is,  two  closely  spaced  monopoles  of  equal  strength  operating  180°  out  of 
phase  constitute  a  dipole  and  two  closely  spaced  dipoles  of  equal  strength 
operating  180°  out  of  phase  with  each  other  constitute  a  quadrupole. 

The  sharp  changes  in  Real  (v‘)  in  the  flame  region  (see  Figures  8  and  9) 
indicate  a  monopole  type  of  acoustic  source.  The  periodic  expansion  of  the 
gases  as  they  move  through  the  flame  under  the  influence  of  an  acoustic  field 
is  similar  to,  in  notion,  to  that  of  a  periodically  expanding  and  contracting 
balloon.  If  this  periodic  expansion  and  contraction  of  the  gases  is  in  phase 
with  the  pressure  oscillations  then  the  resulting  pumping  action  feeds  energy 
into  the  acoustic  oscillations.  If  the  periodic  expansion  and  contraction  is 
out  of  phase  with  the  acoustic  oscillations  energy  is  lost  from  the  acoustic 
motions. 

However,  from  considerations  of  the  spatial  distribution  of  the  heat 
release  rate  fluctuations  it  was  found  that  this  picture  was  not  complete. 
Experimentally,  the  heat  release  rate  fluctuations  may  be  obtained  in  terms 
of  emitted  CH  or  CC  radiation  from  the  flame  (see  Ref.  9).  These  radiation 
emissions  were  measured  as  a  function  of  y  by  capturing  the  radiation  from 
narrow(2  mm  wide)  slits  aligned  perpendicular  to  the  normal  coordinate  y  by 
means  of  a  photomultiplier  arrangement.  As  expected,  close  to  the  sidewall, 
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Developed  Experimental  Set- 


MEASURED  MODEL  INPUTS 


the  measurements  captured  only  the  shot  noise  from  the  photomultiplier 
indicating  no  ‘  gnificant  heat  release  rates  in  this  region. 

The  interesting  discoveries  were  made  in  the  vicinity  of  the  steady 
flame  location.  Just  upstream  of  this  location,  a  large  peak  in  the  magnitude 
of  the  unsteady  heat  release  rate  was  obtained.  At  the  location  of  the  steady 
flame  hardly  any  unsteady  radiation  could  be  detected  which  at  first  glance 
may  seem  surprising  but  is  not  so  as  will  be  explained  shortly.  Just  down¬ 
stream  of  this  location,  however,  the  radiation  levels  peaked  again  although 
the  peak  was  smaller  than  the  first  observed  peak  just  upstream  of  the  steady 
flame  location.  This  behavior  was  found  to  occur  at  frequencies  up  to  1000 
Hz.  When  the  theoretical  model  was  applied  to  obtain  the  magnitude  of  the 
unsteady  heat  release  rates  as  a  function  of  y,  it  showed  identical  results. 
An  example  is  shown  in  Fig.  12,  which  plots  the  theoretical  radiation  levels 
(normalized)  as  a  function  of  y  for  an  arbitrary  driving  frequency.  The 
corresponding  experimental  result  is  shown  in  Fig.  13. 

In  addition,  the  model  predicted  that  for  all  frequencies  the  heat 
release  rate  fluctuations  upstream  and  downstream  of  the  steady  flame  loca¬ 
tion  differed  in  phase  by  180°  (Fig.  14).  According  to  Rayleigh's  criterion 
the  heat  release  rates  may  be  considered  as  the  acoustic  source  of  the  flame. 
These  peaks  on  either  side  of  the  steady  flame  location  are,  according  to  the 
model,  180°  out  of  phase  with  each  other.  If  the  magnitude  of  the  upstream 
peak  is  thought  of  as  the  sum  of  two  parts,  one  of  magnitude  equal  to  the 
downstream  peak  and  the  other  accounting  for  the  remainder,  then  theoretical¬ 
ly  the  flame  may  be  viewed  as  a  combination  of  an  acoustic  dipole  and  a 
monopole.  This  result  is  entirely  new  and  to  the  best  of  the  investigators 
knowledge,  not  found  in  the  literature.  This  split  is  important  from  the 
point  of  view  of  the  directionality  of  the  sound  emitted  by  the  flame.  While 
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Fig.  14.  Predicted  Values  of  the  Phase  of  the  Unsteady  Heat  i 
Rate  (Radiation)  with  Respect  to  the  Imposed  Pressu 
Oscilla  tion. 


MEASURED  OSCILLATORY  RADIATION 

(normalized  magnitudes  and  phases) 


Fig.  13.  Measured  Oscillatory  Radiation  (Heat  Release  Rate)  in  the 
Flame  Region. 
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Fig.  14.  Predicted  Values  of  the  Phase  of  the  Unsteady  Heat  Release 
Rate  (Radiation)  with  Respect  to  the  Imposed  Pressure 
Oscilla  tion. 


a  monopole  source  is  non  directional  in  nature,  a  dipole  source  shows  a 
strong  preference  for  channeling  the  emitted  sound  waves  along  its  ax- 
is(determined  by  joining  the  centers  of  the  two  monopoles  constituting  the 
dipole).  Thus,  for  a  ducted  flame,  the  monopole  source  will  exhibit  a  direc¬ 
tional  preference  purely  due  to  wave  guide  effects  created  by  the  confining 
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effects  of  the  duct  walls  while  the  dipole  source  introduces  additional 
directionality  of  its  own.  This  may  become  important  if  transverse  acoustic 
modes  are  considered. 

The  only  discrepancy  between  the  theoretical  model  and  the  experimental 
measurements  was  found  in  the  phase  difference  between  the  radiation  (heat 
release  rate  fluctuation)  peaks  upstream  and  downstream  of  the  flame.  Experi¬ 
mentally,  it  was  found  that  the  phase  difference  was  not  always  180°  but 
varied  between  130°  to  180.°  This  may  not  be,  however,  entirely  the  model's 
fault  as  the  spatial  resolution  of  the  radiation  measurements  was  limited  by 
data  acquisition  constraints.  A  set  of  experimentally  obtained  phase  values 
is  also  presented  in  Fig.  13. 

The  reason  why  the  radiation  fluctuations  at  the  steady  state  flame 
location  are  small  in  comparison  to  the  sharp  peaks  obtained  upstream  and 
downstream  (Figures  12  and  13)  remains  to  be  discussed.  This  may  be  under¬ 
stood  by  considering  Fig.  11.  At  any  location  y,  the  magnitude  of  the  radia¬ 
tion  fluctuations  is  equal  to  the  difference  in  the  values  of  the  levels  at  t 
=  0  and  t  =  0.5T.  Upstream  and  downstream  of  the  steady  flame  location 
(situated  in  the  flame  displacement  zone)  this  difference  is  seen  to  be  high 
with  the  upstream  levels  being  of  greater  magnitude.  This  corresponds  to  the 
two  peaks  noted  earlier.  However,  in  the  region  close  to  the  steady  state 
flame  location,  the  curves  of  the  instantaneous  heat  release  rates  at  all 
times  during  a  cycle  of  oscillation  intersect.  This  indicates  very  weak 
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fluctuations  in  the  reaction  rate  at  this  location  and  explains  the  low 
radiation  fluctuation  levels  obtained  here. 

Summary 

Keeping  in  mind  the  previously  stated  objectives  of  Task  I  of  the 
research  program  the  following  conclusions  may  be  drawn: 

(i)  The  close  agreement  between  theory  and  experiment  indicates  that 

state  of  the  art  models  of  unsteady  solid  propellant  flames  are 
indeed  capable  of  predicting  the  flame  acoustic  interactions  during 
axial  instabilities. 

(ii)  The  effect  of  the  flame  driving  or  damping  characteristics  is 

manifested  in  sharp  changes  in  the  normal  velocity  fluctuations  in 
the  flame  region.  These  normal  velocity  oscillations  are  equivalent 
to  a  piston  which  supplies  or  receives  acoustic  energy  from  the 
core  flow  oscillations. 

(iii)  Unsteady  heat  transfer  from  the  flame  to  the  propellant  surface  may 
provide  a  feedback  mechanism  between  the  propellant  burn  rate  and 
acoustic  oscillations  which  supplies  the  energy  needed  to  sustain 
axial  instabilities. 

(iv)  The  acoustic  admittance  at  the  propellant  surface  may  exert 
considerable 

influence  upon  the  overall  driving/damping  characteristics  of  the 
rocket  motor. 

(v)  The  flame  may  be  considered  as  a  combination  of  an  acoustic  monopole  and 
an  acoustic  dipole.  The  dipole  nature  may  become  important  in  studies  of 
transverse  wave  mode  instabilities 
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Task  II 


Interaction  of  Side  Wail  Stabilized  Diffusion 
Flames  with  Axial  Acoustic  Fields 

Task  II  of  this  research  program  investigated  the  role  of 
diffusion  limited  processes  upon  the  driving/damping  of  axial  instabilities 
in  solid  propellant  rocket  motors.  The  major  objectives  of  this  task  were: 

1.  To  incorporate  the  effect  of  diffusion  processes  in  theoretical  models 
unsteady  solid  propellant  flames  and 

2.  To  investigate  the  influence  of  diffusion  processes  in  the  flame  region 
on  the  acoustic  driving/damping  characteristics  of  the  flame. 

The  utilized  experimental  setup  is  described  first.  This  is  followed  by 
a  discussion  of  the  conducted  theoretical  and  experimental  investigations. 

Experimental  Setup 

A  schematic  of  the  developed  experimental  set  up  is  shown  in  Fig.  15. 

It  consists  of  a  long  rectangular  tube  having  a  diffusion  flame  burner  on  one 
of  its  side  walls,  an  axially  movable  injector  plate  at  the  inlet  end  and  two 
phase  locked  acoustic  drivers  at  the  exhaust  end.  During  an  experiment, 
diffusion  flames  are  stabilized  on  the  lower  burner  surface.  Next,  the 
acoustic  drivers  are  turned  on  to  excite  a  standing  longitudinal  acoustic 
wave  of  desired  frequency  and  amplitude  in  the  duct.  The  latter  simulates 
the  axial  oscillations  observed  in  unstable  solid  propellant  rocket  motors. 
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Schematic  of  the  Proposed  Experimental  Set-Up  to  Investigate  the 
Interactions  between  Gas  Phase  Diffusion  Flames  and  Longitudinal 
Acoustic  Fields. 


The  position  of  the  burner  and,  consequently  the  positions  of  the  flames 
relative  to  the  standing  acoustic  wave  (i.e.,  next  to  a  pressure  node  or 
pressure  antinode)  can  be  varied  by  moving  the  injector  plate  axially.' 

The  diffusion  flame  burner  utilized  during  the  studies  conducted  during 
the  past  year  is  shown  in  Fig.  16.  Fuel  (propane)  and  air  are  injected 
through  alternate  slots  and  this  results  in  three  diffusion  flames  being 
established  at  the  burner  surface.  This  type  of  burner  has  two  advantages. 
First,  the  resulting  diffusion  flame  configuration  is  geometrically  simple 
and,  therefore,  amenable  to  theoretical  analysis.  Secondly,  the  developed 
diffusion  flames  are  similar  to  the  sandwich  type  of  propellant  flames  with 
alternating  oxidizer  and  binder  sections  (see  Fig.  2)  which  have  been  uti- 
lized  in  the  past  to  study  solid  propellant  combustion  . 

Theoretical  Investigations 

One  of  the  objectives  of  the  diffusion  flame  studies  is  the  development 
of  a  theoretical  model  of  oscillatory  diffusion  flames  stabilized  in  longitu¬ 
dinal  sound  fields  simulating  those  found  in  unstable  solid  propellant  rocket 
motors.  Such  a  model,  based  upon  the  Schvab-Zeldovich  coupling  function 
formalism1*,  has  been  developed  during  this  program.  This  approach  was 
adopted  as  it  is  known  to  be  capable  of  predicting  the  behavior  of  diffusion 
flames  under  steady  conditions.  Moreover,  it  captures  the  effects  of  diffu¬ 
sion  processes  which  are  the  focus  of  this  phase  of  the  investigation. 

The  Schvab-Zeldovich  approach  is  most  applicable  to  those  cases  in  which 
the  diffusion  rates  of  fuel  and  oxidizer  towards  each  other  control  the 
overall  reaction  rates;  that  is,  cases  where  the  chemical  reaction  tim^  is 
very  much  smaller  that  the  time  required  for  the  diffusion  of  the  reactants 


18 


0  - 1  cm 


re  Gauge 


Ceramic  Matrix 


Ceramic  Matrix 


Fuel  Manifold 


Sintered  Steel  Pipe 


Air 


Fig  16. 


$ 


loped  Diffusion  Flame  Burner. 


Configuration  for  the  Theoretical  Model  of  the  Diffusion  Flame 
Acoustic  Interactions. 


towards  each  other.  In  such  cases,  the  actual  flame  region  may  be  taken  to 
lie  on  a  surface  or  sheet.  Moreover,  by  appropriate  definitions  of  tempera¬ 
ture  and  species  coupling  functions, ft  ^  and  fi ,  respectively,  the  chemical 
reaction  rates  may  be  eliminated  from  consideration. 

The  theoretical  model  considers  the  configuration  shown  in  Fig.  17.  A 
number  of  diffusion  flames  are  established  on  the  side  wall  of  a  duct  by 
means  of  alternate  fuel  and  oxidizer  slots.  This  configuration  is  similar  to 
the.  experimental  set  up  considered  in  Figs.  15  and  16.  These  flames  are 
exposed  to  known  axial  acoustic  pressure  and  velocity  oscillations.  Of 
interest  is  the  response  of  these  flames  to  the  imposed  oscillations.  It  is 
sufficient  to  consider  a  single  flame  as  the  response  of  any  one  of  the 
flames  is  largely  governed  by  the  conditions  existing  locally.  Thus,  the 
model  considers  the  single  diffusion  flame  shown  in  Fig.  17.  The  influence 
of  the  other  flames  is  accommodated  by  prescribing  the  flow  variables  along 
the  boundary  of  the  domain  of  interest  as  shown  in  Fig.  17. 

The  model  considers  three  unknowns;  namely,  (1)  the  normal  velocity 
fluctuation  v‘  (2)  the  fluctuations  in  the  temperature  coupling  function  U'j, 
and  (3)  the  fluctuations  in  the  species  coupling  function,  ' . 

These  unknowns  are  related  to  the  applied  acoustic  field  described  in 
terms  of  known  pressure  fluctuations,  p1,  and  axial  velocity  fluctuations, 
u1,  by  making  use  of  the  conservation  equations  for  mass,  momentun  and  energy 
along  with  the  equation  of  state.  A  matrix  equation  of  the  following  form 
emerges : 

X,+  C,X,=  F: 

where 

X=(v\  pT,  pV 
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and  A..  B .  .  and  C • •  are  3x3  coefficient  matrices  which  depend  upon  the 

'  J  5  '  J  '  J 

steady  state  solutions.  These  have  to  be  obtained  experimentally  or  by  a 
combination  of  theory  and  experiment.  F.  depends  upon  p'  and  u'  as  well  as 
upon  the  steady  state  solutions  and  may  be  regarded  as  a  forcing  function 
for  the  unknowns  X..  Conditions  on  the  velocity,  temperature  and  species 
fractions  (i.e.,  on  v‘,/3'^.  and  /3' )  have  to  be  prescribed  on  the  boundary  as 
depicted  in  Fig.  7.  These  boundary  conditions  have  to  be  obtained  experimen¬ 
tally  for  a  general  case.  However,  these  conditions  may  be  specified  approx¬ 
imately  as  follows: 

(a)  At  the  burner  surface  (i.e,  y=0)  the  fuel  and  oxidizer  streams  are 
each  honmogeneous  and  separate  from  one  another  so  that  there  are  no  fluctua¬ 
tions  in  the  species  fractions  at  this  location.  Moreover,  temperature 
fluctuations  at  the  ceramic  matrix  surface  of  the  burner  are  expected  to  be 
negligible.  This  implies  that  at  y=0,  &' j  =0  and  P  -  0.  The  condition  on 
v1  at  y=0  is  specified  in  terms  of  the  acoustic  admittance  R  of  the  burner 
surface;  that  is, 

v'(y=0)  =  Rp* 

(b)  At  x=0  and  at  x=l,  conditions  may  be  expected  to  be  symmetric  if  one 
considers  the  middle  diffusion  flame  (which  has  similar  diffusion  flames  on 
either  side.)  Hence,  the  conditions  may  be  specified  to  be 

dPT  /  dx  =  Q  and  dpi  dx=0  at  x=0  and  at  x  =1 

A  computer  program  for  solving  the  model  equations  has  been  developed. 
A  Crank-Nicholson  implicit  numerical  scheme  which  employs  uniform  grid 
spacings  in  both  the  x  and  y  directions  is  utilized  to  obtain  a  finite 
difference  representation  of  the  developed  matrix  equation.  A  block 
tridiagonal  system  of  finite  difference  equations  is  formed  at  each 
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y-location  which  can  be  inverted  efficiently  using  standard  algorithms  to 
obtain  the  required  unsteady  solutions. 

Results  obtained  using  a  hypothetical  set  of  input  and  forcing  condi¬ 
tions  are  shown  in  Figs.  18  and  19.  For  this  case,  the  burner  surface 
admittance  has  been  taken  to  be  purely  reactive;  that  is,  the  real  part  of  v1 
(the  component  in  phase  with  p‘)  at  the  burner  surface  is  zero  and  v'  leads 
p'  by  90°  at  the  surface.  Figure  18  describes  the  variation  of  the  magnitude 
of  v'  at  the  top  edge  of  the  flame  (i.e.,  at  y  =  y^)  as  a  function  of  x  (see 
Fig.  17 )>  Figure  19  describes  the  variations  of  the  phase  of  v*  with  respect 
to  p'  at  this  location  (i.e.,  y=y^)  as  a  function  of  x.  These  results  show 
that  for  the  assumed  conditions  the  magnitude  of  v'  exhibits  a  maximum  in  the 
vicinity  of  the  flame  surface.  Figure  19  also  shows  that  the  phase  differ¬ 
ence  between  p'  and  v‘  is  less  than  90°  in  this  region.  This  indicates  the 
presence  of  a  positive  real  part  of  v'  (i.e.,  the  component  in  phase  with 
p').  As  shown  previously,  the  product  of  the  in-phase  components  of  p'  and 
v'  yields  the  energy  input  into  the  acoustic  wave  by  the  flame  in  the  vicini¬ 
ty  of  the  flame  surface.  A  smaller  peak  in  the  magnitude  of  v'  appears  a 
short  distance  axially  Jownstream  of  the  flame  surface.  Significantly,  in 
this  region  (see  Fig.  19),  p1  and  v'  are  out  of  phase  indicating  a  damping  of 
the  acoustic  motions  and  the  presence  of  a  negative  real  part  of  v1  in  this 
region. 

Thus,  the  model  predicts  the  occurrence  of  driving  and  damping  regions 
within  the  flame.  This  is  similar  to  the  monopole-dipole  mechanism  occurring 
in  the  premixed  flame  case  which  was  investigated  during  Phase  I  of  this 
program.  This  similarity  suggests  that  the  developed  diffusion  flame  model 
may  indeed  be  capable  of  predicting  the  driving  characteristics  of  diffusion 
flames. 
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Experimental  Efforts 

The  experimental  studies  conducted  to  date  included  high  speed  flame 
shadowgraphy  visualizations  and  measurements  of  CH  radiation  from  the  flame. 

High  speed  shadowgraph  films  (5000  frames/sec)  have  been  taken  with  the 
flame  excited  at  different  frequencies  in  the  250  -  1000  Hz.  range  by  means 
of  the  acoustic  drivers  (see  Fig.  15).  The  axially  movable  end  plate  made  it 
possible  to  place  the  diffusion  flame  burner  on  different  locations  of  the 
established  standing  wave;  that  is,  at  a  pressure  maximum,  minimum  or  in 
between.  Analyses  of  the  films  indicate  that  under  the  influence  of  an 
acoustic  field  the  flame  oscillates  axially  and  the  frequency  of  oscillation 
coincides  with  the  frequency  of  the  excited  acoustic  wave. 

Measurements  of  C-H  radiation  emitted  by  the  diffusion  flames  have  also 
been  carried  out  for  different  frequencies  and  flame  locations  as  noted 
above.  These  measurements  are  indicative  of  the  heat  release  rates  of  the 
flames.  An  example  of  the  measured  phase  difference  between  the  flame 
radiation  and  acoustic  pressure  oscillations  at  a  pressure  maximum  for 
different  frequencies  in  the  250  to  1000  Hz.  range  is  shown  in  Fig.  20. 
According  to  Rayleigh's  criterion,  when  the  phase  difference  between  the  two 
is  less  than  90°  the  flame  drives  the  acoustic  field.  Thus,  it  may  be  noted 
that  depending  upon  the  frequency,  the  flame  may  either  drive  or  damp  the 
acoustic  oscillations. 

Due  to  time  constraints,  the  validity  of  the  developed  theoretical  model 
could  not  be  completely  assessed  during  the  period  of  this  grant.  However, 
its  validity  will  be  ascertained  under  an  extension  of  this  program  in  which 
the  major  focus  will  shift  to  the  importance  of  flow  turning  phenomena  in 
damping  and  driving  axial  instabilities.  The  model  will  be  checked  by 
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comparing  the  distribution  of  the  unsteady  velocity  v'  measured  using  an  LDV 
with  the  theoretical  predictions.  Measurements  of  the  steady  state  tempera¬ 
ture  distributions  in  the  flame  region  and  the  acoustic  admittance  of  the 
ceramic  burner  surface  will  also  be  needed  as  inputs  to  the  theoretical 
model . 


Summary 

The  studies  conducted  under  Task  II  of  the  program  have  shown: 

(i)  Side  wall  stabilized  diffusion  flames  are  responsive  to  axial 
acoustic  oscillations 

(ii)  The  acoustic  driving/damping  characteristics  of  diffusion  flames 
are  frequency  dependent. 

(iii)  The  acoustic  nature  of  diffusion  flames  is  similar  to  the 
monopole-dipole  characteristics  of  premixed  flames. 

(iv)  Diffusion  processes  may  be  incorporated  into  theoretical  models 
which  describe  unsteady  solid  propellant  flame  behavior. 
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